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ABSTRACT 

We detected three extremely embedded X-ray sources in the R Corona Aus- 
tralis (R CrA) star for min g core, IRS 7 region. Two weak X-ray sources are 
associated with the VLA centimeter radio sources 10E & W, whereas the third 
brightest source detected in the two XMM-Newton observations on March 2003 
has no counterpart at any wavelengths. The large R-band upper- limit (19.4 m ) 
measured with the University of Hawaii 88-inch Telescope and strong absorption 
derived in X-rays (Ah ~2.8x10 23 cm -2 equivalent to A\ ~180 m ) indicate that 
the source is younger than typical Class I protostars, i.e. a Class 0 protostar or 
an intermittent phase between Class 0 and Class I protostars. The X-ray lumi- 
nosity was less than one thirtieth (log Lx <29.3 ergs s -1 ) in the former Chandra 
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observation in October 2000, which suggests that the X-ray activity, probably 
generated by magnetic activity, is triggered by an intermittent mass accretion 
episode such as FU Ori type outbursts. Because the source was detected at high 
significance in the XMM-Newton observations (~2,000 cnts), X-ray properties 
of such young protostars can be well investigated for the first time. The light 
curves were constant in the 1st observation and increased linearly by a factor of 
two during 30 ksec in the 2nd observation. Both spectra showed iron K lines 
originated in hot thin-thermal plasma and fluorescence by cold gas. They can be 
reproduced by an absorbed thin-thermal plasma model with a Gaussian compo- 
nent at 6.4 keV (kT ~3— 4 keV, Z-x ~7— 20xl0 30 ergs s _1 ). The rising timescale 
of the light curves in the 2nd observation was too slow for magnetically generated 
X-ray flares, whereas large equivalent width of the fluorescence iron K line in the 
1st observation (~810 eV) requires strong partial covering of the X-ray source. 
These results suggest that a confined hot (perhaps accretion) spot on the proto- 
stellar core was behind the star in the 1st observation and just appeared in the 
2nd observation due to the core rotation with period of >2.8 days, which is much 
slower than the break-up velocity previously assumed for young protostars. This 
means that the source had quiescent X-ray activity during the observations with 
an order of magnitude stronger level than Class I (older) protostars. We also 
consider whether the X-ray source associated with 10E could be shock heated 
plasma by a collision of a jet emanating from the tentative Class 0 object. 

Subject headings: stars: activity — stars: pre-main-sequence — X-rays: stars 


1. Introduction 

The youngest phase of solar type stars, low mass protostars, are found deep inside molec- 

ular cloud cores. Violent mass accretion is believed to occur onto the proto-stellar core, but 

the activity is generally hidden inside an enormous thickness of gas. From spectral energy 

distribution (SED) in the infrared (IR) and radio wavelengths, protostars are divided into 
two classes, which are thought to represent their evolutionary status. The former phase of 
protostars designated as Class 0 objects (~10 1 * * 4 -year-old) mainly have far-IR and submillime- 
ter emission from cold circumstellar dust (T <30 K, Andre et al. 1993). They are thought 
to undergo dynamical mass accretion from huge circumstellar envelopes. The latter phase 
of protostars designated as Class I objects (~10 5 -year-old) can be observed in the near-IR 
wavelength. The proto-stellar cores become hot (T ~3000— 5000 K) by assembling most of 
their stellar mass (late accretion phase). 
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X-ray observations, especially in the hard X-ray band transparent for thick molecular 
clouds, have shown high energy activity on Class I objects (Koyama et al. 1996; Grosso 
et al. 1997; Imanishi et al. 2001). The X-ray activity is probably driven by magnetic activity 
close to the proto-stellar core as is suggested by occasional rapid X-ray flares reminiscent of 
solar flares. A solar-type magnetic dynamo does not work on protostars according to current 
theories. As an alternative mechanism, reconnection of fossil magnetic fields linked to the 
mass accretion process has been proposed (Koyama et al. 1996; Montmerle et al. 2000). 

X-ray emission from extremely embedded sources (Ah ~1— 3xl0 23 cm -2 ) was first re- 
ported in the OMC-2/3 cloud (Tsuboi et al. 2001). They have characteristics of Class 0 object 
as having no near-IR counterparts and associating with millimeter radio clumps. However, 
follow-up radio and near-IR observations by Tsujimoto et al. (2004) did not classify them 
as Class 0 protostars unambiguously. In particular, one of the two candidates have good 
correlation with a centimeter radio source and a jet feature in the H 2 band, evidence for 
excitation by a jet ejected from a nearby Class I protostar. Furthermore, X-ray photons ob- 
tained in the observation (~50 photons) was too small to deduce reliable plasma parameters 
or variability. Though there is also other X-ray detections from cores associated with Class 0 
candidates in the Trifid nebula (Rho et al. 2004), long distance to the sources (d 1.7 kpc) 
makes their identification fax more difficult. To date no X-ray sources has been identified 
with a bona-fide Class 0 object. 

The R Corona Australis (R CrA) cloud is a nearby low and intermediate-mass star 
forming region ( d ~170 pc, Knude &; H<pg 1998). Among many young stellar objects in the 
cloud, those in the IRS 7 region have attracted particular interest as a cite of ongoing star 
formation. The region contains double peaked strong centimeter emission; one in the east 
designated as 10E (or IRS 7B) and the other in the west designated as 10W (IRS 7A) (Brown 
1987; Feigelson et al. 1998; Harju et al. 2001). It also has corresponding two submillimeter 
peaks (van den Ancker 1999), multiple millimeter continuum peaks (Henning et al. 1994; 
Saraceno et al. 1996; Chini et al. 2003; Choi & Tatematsu 2004) and signature of strong 
bipolar outflows (Harju et al. 1993; Anderson et al. 1997), whereas it has only one near-IR 
source associated with 10W (Wilking et al. 1997). These characteristics make the IRS 7 
region a promising host of Class 0 sources. Koyama et al. (1996) detected hard X-ray 
emission and an intense flare from the IRS 7 region, possibly from IRS 7 (Rl), suggesting a 
protostar with high energy (magnetic) activity. 

In this paper, we report on the X-ray detection of extremely embedded sources in the 
IRS 7 star forming core dining XMM-Newton observations in 2003. These results combined 
with analyses of an earlier Chandra observation in 2001 and follow-up near-IR observations 
with the University of Hawaii 88-inch Telescope (UH88) help to understand the nature of 
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the detected X-ray sources. 


2. X-ray Observations and Absolute Position Correction 

We analyzed the X-ray data obtained with the XMM-Newton observatory in two ob- 
servations on 2003 March 28 (hereafter ObsxMMi) and March 29 (ObsxMM2) (Table 1). The 
XMM-Newton satellite (Jansen et al. 2001) has three nested Wolter I type X-ray telescopes 
(Aschenbach et al. 2000) with two types of CCD cameras (one EPIC-pn camera and two 
EPIC MOS cameras, Striider et al. 2001; Turner et al. 2001) on their focal planes. During 
the observations, all CCD cameras were operated in full frame mode with the medium filter 
inserted. Obs X MMi and ObsxMM2 had the same satellite pointings. The IRS 7 star forming 
core was at 6' off-axis, where the 90% photon radius is ~1'. The analysis was performed 
with the software package SAS ver. 5.4.1 and LHEASOFT ver. 5.2. The observing data 
were processed using the SAS scripts “emchain” and “epchain” . Events detected during pn 
frame transfer, called “out of time events”, were reproduced as well. Events with pattern 
^12 (pattern ^4 for pn spectral analysis) were selected. High background periods, defined 
as count rate of “pattern = 0” events above 10 keV across the entire chip above 0.35 cnts s -1 
for MOS and 1 cnts s -1 for pn, were screened out. Events close to hot pixels or outside the 
field of view (FLAG = 0) were also removed. 

To help source identification and investigate long-term flux variations, we analyzed the 
archival Chandra observation on 2000 October 7 (hereafter Obscxo, Table 1), whose first 
result was presented by Garmire & Garmire (2002). Chandra has one X-ray telescope with 
capability of sub-arcsecond imaging and two inter-changeable detectors (Weisskopf et al. 
2002). In the observation, the ACIS-I CCD camera arrays were selected. The IRS 7 star 
forming region was at 2' off-axis, where the 90% photon radius is ~1." 3. The analysis was 
performed with the software package CIAO ver. 2.3 and LHEASOFT ver. 5.2. The level 1 
event data were reprocessed with “acis_process_events” to compensate for degradation of the 
CCD spectral resolution by charge transfer inefficiency (CTI, Townsley et al. 2000). Events 
with grade 0, 2, 3, 4 and 6 were selected. 

For correction of absolute coordinates, positions of nine X-ray bright sources (two 
more sources for the Chandra data) on the EPIC pn plus MOS image (0.3—10 keV) and 
the ACIS-I image (0.5—9 keV) were measured using source detection packages (SAS: ede- 
tecLchain, CIAO: wavdetect). These positions are correlated with 2/xm IR sources in the 
2MASS All-Sky Point Source Catalog (PSC) obtained from NASA/IPAC Infrared Science 
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Archive (GATOR) 1 . The astrometric accuracy of this catalogue is good around 0." I 2 . Each 
X-ray source has only one 2MASS counterpart within 4", so that source misidentification 
is not a problem. After adjusting the X-ray data frames as shown in Table 1, deviation of 
Chandra and XMM-Newton coordinates from the 2MASS frame is < 0." 6. 


3. X-ray Image 

During the XMM-Newton observations, strong hard X-ray emission was detected from 
the IRS 7 star forming core. Figure 1 shows the true color EPIC pn + MOS image of 
the R CrA star forming region combining the ObsxMMi and ObsxMM2 data. The image is 
color-coded to represent soft band (0.2—1 keV) to red, medium band (1—3 keV) to green 
an d hard band (3—9 keV) to blue. The IRS 7 region has darker blue than the other Class I 
candidates, suggesting that the emission is heavily absorbed. The left panel of Figure 2 
magnifies the IRS 7 region of the blue band (3—9 keV) image, shown in black and white. A 
strong X-ray source (X| MM ) was detected near the radio source 10E in both ObsxMMi and 
Obs X MM 2 (Table 2). The coordinates derived from each observation with edetect-chain are 
consistent (<0." 3 discrepancy) so that the sources seen in ObsxMMi and ObsxMM 2 would 
be identical. We also found another faint X-ray source (X™ M , ~14a significance for the 
ObsxMMi data) to the north-east direction of Xg MM . The source was not detected with 
edetecCchain probably due to crowdedness of the region. 

Chandra detected two weak X-ray sources in the IRS 7 region at above 6(7 significance 
(Xg xo and X& xo , see Table 2 and the top right panel of Figure 2). X^ xo is identified 
with X™ M (hereafter we call it Xw), but Xg XO has a 2" offset to south-east from Xg MM 
although position accuracy of Xg MM could have additional unknown uncertainty affected by 
unresolved sources. This suggests no X-ray emission from Xg MM detected in Obscxo? while 
Xg XO emission in ObsxMMi and Obs X MM2, if present, can be obscured by the wings of X| MM 
point-spread-function. (We think that the component of Xg XO is marginally resolved in the 
spectra of Xg MM , see section 4.1). 

Compared with field sources in the other wavelengths (the bottom right panel of Fig- 
rue 2), X w has VLA centimeter (10W in Brown 1987), millimeter (source 4 in Choi & 
Tatematsu 2004) radio, and neax-IR (IRS 7) counterparts. Xg XO has a VLA centimeter 
(10E in Brown 1987) radio counterpart. Xf MM has no counterpart in any other wavelength. 


1 http://irsa.ipac.caltech.edu /applications/Gator/ 

2 Sect. 2.2 of Explanatory Supplement to the 2MASS All-Sky Data Release: 

http://www.ipac.caltech.edu/2mass/releases/allsky/doc/ 
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4. X-ray Light Curves and Spectra 

With ~2,000 X-ray events obtained from Xg MM with XMM-Newton, we could perform 
detailed timing and spectral analysis. Photon statistics events of the other sources are quite 
limited (10—120 photons). We briefly show their spectral analyses. 


4.1. East Sources (Xg MM & Xg xo ) 

Source events of X^ MM were extracted from a 27." 5 radius circle centered on the source 
position excluding the region around X((X M (like a crescent shape). Background was selected 
from source free regions. Source events from X£ XO were extracted from a 1." 8 radius circle 
centered on the source position. The background has a negligible contamination. 

Plots in the top panel of Figure 3 show background subtracted light curves of X^ MM in 
the 2—10 keV band. The light curves in ObsxMMi (first half) were rather flat with marginal 
fluctuation in the latter observation. A constant model in a x 2 fit is acceptable at the 90% 
confidence (Table 3). Count rates in ObsxMM 2 started at around four times as high as the 
average count rate in ObsxMMi • The light curves gradually increased by a factor of two 
throughout ObsxMM2 except for two dips at around 8.02 and 8.17 days. They are consistent 
with a linear model with a slope of ~9.3 xl0~ 2 cnts s -1 day -1 . The hardness ratio defined 
as count rates in the 5—10 keV band over those in the 2—5 keV band (bottom panel of 
Figure 3) remained unchanged within statistics both in ObsxMMi and ObsxMM 2 except for 
a strange minimum at 7.86 day. This means that time variation is not accompanied by 
significant change of the spectral shape. On the other hand, the X-ray dips on timescales of 
<1 ksec are unusual, but they do not seem to be artificial: both pn and MOS light curves 
without background subtraction also show such dips; other nearby sources (R CrA, IRS 5 
and CrA 1) have no dips at those timings; similar dips are seen at 7.85 (only MOS data 
shown with diamonds in Figure 3 is available) and 6.56 days with less significance (see arrows 
in Figure 3). Interestingly, these dips are on a period of ~13.9 ksec with the exception of 
the dip at 6.40 days which is missing. 

The spectra of X| MM in ObsxMMi and ObsxMM 2 (Figure 4) show several similarities. 
X-ray emission in both reaches up to ~10 keV and is strongly absorbed below 2—3 keV. 
It has a broad line feature between 6—7 keV, and marginal lines, which are not obvious in 
the EPIC MOS spectra, between 5—6 keV. Origin of the marginal line feature is a puzzle 
because there is no emission line of any abundant element between 5—6 keV, whereas the 
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background chromium line of the EPIC pn at 5.5 keV 3 is too weak to show such a feature. 
To investigate the line feature between 6-7 keV, we tentatively tried a simultaneous fit of the 
pn and MOS spectra of ObsxMM 2 by an absorbed power-law model with a Gaussian compo- 
nent. An acceptable fit at above 90% confidence has the photon index of 3.0 (2. 5-3. 4), the 
Gaussian center energy of 6.60 (6.53-6.67) keV and the Gaussian sigma of 0.15 (0.079-0.28) 
keV (braces show 90% confidence range). Taking into account the source as a promising 
young stellar origin (see discussion 6.1), the significant line broadening is reasonably re- 
garded as a blend of line emission from He- like iron in hot plasma (~6.7 keV) and neutral 
(fluorescent) iron (~6.4 keV) in cold gas rather than a Doppler broadening with Av ~7,000 
km s -1 . We therefore simultaneously fit both pn and MOS spectra for each observation 
by an absorbed 1-temperature (IT) thin-thermal plasma model (absorption: Morrison & 
McCammon 1983, emission: the MeKaL plasma code, Mewe et al. 1995) with a Gaussian 
component with the fixed line center at 6.4 keV. Fits for each observation are acceptable 
above 90% confidence (Table 4, Model A and B for ObsxMMi and Obs XM M2, respectively), 
but A H has significant discrepancy between ObsxMMi and Obs XM M2- Ah for ObsjtMMi is 
determined from the slope between 1—5 keV, whereas that for Obs XMM 2 is from the steep 
slope around 3-5 keV, neglecting weak emission below 3 keV. Considering that emission 
below 3 keV seems unchanged between Obs X MMi and Obs XMM2 in Figure 4, we believe that 
a model with soft constant emission below 3 keV and hard variable emission above is more 
realistic. We thus refit the spectra of Obs XM Mi and Obs XM M2 simultaneously by an absorbed 
2T model with a Gaussian component at 6.4 keV. In the model, a IT component for the soft 
constant emission, A H of the hard variable component and elemental abundance for both 
components were tied together between ObsxMMi and Obs X MM 2 - Whereas, Ah parameters 
of the constant and variable components were independently varied because a model with a 
common Ah gives unrealistically large absorption corrected X-ray luminosity for the const an t 
component as a stellar origin (log L x ~ 35 ergs s -1 ). The model (Model C in Table 4), again, 
successfully reproduced the spectra at above 90% confidence. The absorption column in the 
best-fit model is extremely large ~2.8xl0 23 cm -2 , equivalent to Ay ~180 7Tt (Imanishi et al. 
2001). Plasma temper ature is ~2.7 keV in ObsxMMi and rises to ~4.0 keV in Obs X MM2- The 
absorption corrected X-ray luminosity in Obs X MMi (Jog T x ~30.8 ergs s -1 ) is large, which 
further increases in Obs X MM 2 to log L x ~31.2 ergs s -1 . Metal abundance is sub-solar (~0.2 
solar), which is typical of low mass young stars. The fluorescent iron fine has equivalent 
width (EW) of ~810 (240—1400) eV in Obs X MMi> which is one of the largest among young 
stellar origin (e.g. Imanishi et al. 2001). The line intensity did not change significantly in 
Obs X MM 2 , and therefore EW dropped to ~250 (100—400) eV (braces show 90% confidence 


3 Sect. 3. 3. 7.2 in XMM-Newton Users Handbook, 
http: / / xmm.vilspa.esa.es / external /xmm_user_support / documentation /uhb /index, html 
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range). 

The X-ray spectrum of Xe X 0 (left panel of Figure 5) looks rather flat between 1—10 
keV though its statistics is poor (~20 photons). The shape is, roughly speaking, not like 
the hard variable component of Xg MM , which has deep absorption below 3 keV, but like 
the soft constant component. Indeed, Xe XO has 3 cnts below 3 keV, whose 99% confidence 
range is between 0.39 cnts and 10.15 cnts ( N = 3, B — 0 in Table 3 of Kraft et al. 1991). 
If Xg X0 were identical to the hard variable component of Xg MM , the normalized photon 
counts below 3 keV should amount to only ~0.18 cnts. Remembering that Xg XO has a 2" 
offset from X| MM (see section 3) and is included in the source region of Xg MM , it would be 
identical to the minor (=soft constant) component of Xg MM . We modeled the spectrum of 
Xg xo by an absorbed IT model fixing Nn at 4.2xl0 22 cm 2 , the best-fit value of the soft 
constant component, and varying kT. An acceptable model (Model D in Table 4) has hot 
plasma temperature >4 keV and an X-ray flux of one third of the soft constant component 
ofX XMM . 

During Obscxo i no X-ray photon was detected within the error circle of X£ MM (see 
Figure 2), suggesting that the hard variable component of Xg MM was below the detection 
level (<3 photons at the 95% confidence referring to N = 0, B = 0 in Table 2 of Kraft et al. 
1991). Assuming kT ~ 3 keV and jV H ~2.8xl0 23 cm" 2 , the X-ray luminosity is log L x ^29.3 
ergs s -1 , which is ~l/30 of the hard variable component in Obs X MMi- It might be in an 
inactive (dormant) phase during Obscxo- 


4.2. West Source (Xw) 

Source events of X^ MM (X w during ObsxMMi and ObsxMM 2 ) were extracted from an 
ellipse with 15" long and 10." 5 short axes elongated toward the NNE direction, excluding a 
region with strong contamination from X£ MM . The background events are extracted from a 
symmetrical region of the source region with respect to Xg MM . The Obs X MM 2 data were not 
used due to stronger contamination from X XMM . Source events of X^ xo (X w during Obscxo) 
were extracted from a 1." 8 radius circle. The background has a negligible contamination. 

EPIC pn and MOS spectra of X™ M (right panel of Figure 5) were reproduced by an 
absorbed IT model with strong absorption (Ah ~3.4x 10 23 cm -2 ), hot plasma temperature 
(kT >1.6 keV) and log L x ~ 30.5 ergs s -1 (Model E in Table 4). Because photon counts 
of XCX° was quite poor (~13 photons), it is reproduced by the best-fit model of X™ M just 
by changing its normalization (Model F in Table 4). The X-ray luminosity (log L x ~ 30.0 
ergs s -1 ) is around one third of that during ObsxMMi- 
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5. Follow-up Observation in the Near-IR 

Infrared source catalogues currently available have several shortcomings for our study: 
shallow limiting magnitude (e.g. Ku mit ~ 15 m for 2MASS, 16.5 m for Wilking et al. 1997); 
too large pixel scales to resolve objects from nearby nebulosity (e.g. 1" pix -1 for 2MASS and 
0." 75 pix -1 for Wilking et al. 1997); insufficient positional accuracy (e.g., ±1" for Wilking 
et al. 1997). We therefore analyzed two deep K - band images of the field obtained in August 
1998 and August 2003 using UH88 and the near-infrared imager QUIRC. (Full analysis of 
the data covering ~100 / fov will be published elsewhere.) During the observations, the f/10 
secondary mirror was used, yielding the pixel scale of 0." 1886 pix -1 . The spatial resolution 
of our images are 0." 5 and 0." 8 in FWHM for the data of 1998 and 2003, respectively. The 
analysis was made with IRAF 4 . 

Both images were set to the absolute coordinate by referring to the coordinate of R CrA 
in the 2MASS catalogue. Since the positions of R CrA in the UH88 frames were difficult to 
measure because of saturation, we measured the positions of IRS 9 and derived the relation 
between R CrA and IRS 9 from a high resolution image (0." 14) obtained with the SUBARU 
8.2 m telescope with the adaptive optics system (Nedachi et al. in preparation). The fin al 
positional errors of our images are estimated at ~0." 3 at maximum. Both images look the 
same except for their data quality. The top left panel of Figure 6 shows the image taken in 
1998 in better condition. In the image, we see filamentary and mildly extended reflection 
nebulae of R CrA. To highlight fine structures smaller than 10" scales for aiding source 
detection, we median-filtered the image with a 10" x 10" box and subtracted the resulting 
image (the bottom panel of Figure 6) from the original. The result is showm in the top ri gh t, 
panel of Figure 6. We found no significant emission from Xg MM and Xg XO though we see a 
marginal enhancement of flux at the location of Xg X °. To determine their flux upperlimit, 
we flattened the background (mainly extended nebular emission of R CrA) by subtracting 
the fitted third-order polynomial surfaces from the original images, and measured fluctuation 
of a subset of the image containing 20 x 20 pixels at that location of the sources. The 5- 
a detection upperlimit for both sources at 50% flux enclosure radius were ~19.4 m for the 
1998 image and ~19 m for the 2003 image. The absolute coordinate of IRS 7 measured 
with the SUBARU image (Nedachi et al. in preparation) is determined at (a'2000, 62000) — 
(19hlm55.34s, — 36d57’21.69”) within 0." 3 accuracy. With our new coordinate, IRS 7 falls 
in the potitional error circle of X^ xo , as shown in Figure 2. 


4 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 


- 10 - 


6. Discussion 

6.1. The Nature and X-ray Emission Mechanism of Xg MM 

X| MM has no known IR nor radio counterpart. It could be a background AGN serendipi- 
tously discovered. However, AGNs above the observed X-ray flux (2-8 xl0~ 13 ergs cm -2 s -1 
between 2-10 keV) are found <10 degree ~ 2 in the sky (Ueda et al. 1998), and the proba- 
bility to detect such an AGN in the IRS 7 star forming core (— lOarcsec 2 ) is extremely small 
(<10 -4 ). Furthermore, the spectra show thermal iron K emission line unusual for an AGN 
and, if we try to fit the X-ray spectrum in Obs X MM2 by an absorbed power law model, its 
spectral photon index of —3.0 is far larger than the AGN standard (T < 2 for AGNs, Ueda 
et al. 1998). These facts exclude an AGN origin. Taking its location into account, X| MM is 
most likely an unknown very young stellar object. 

The large X-ray absorption column (Ah ~2.8x 10 23 cm -2 ) and small near-IR luminosity 
( K >19. 4 m ) compared with nearby Class I objects such as IRS 1, 2 and 5 (Ah — 2x 10 22 cm -2 : 
our additional analysis of the XMM-Newton data, K <ll m : Wilking et al. 1997) mean that 
X| MM has a huge circumstellar envelope which has not well heated up yet. Certainly, van 
den Ancker (1999) indicated a strong submillimeter (450 &850 /rm) condensation around 
Xg MM . This strongly suggests that Xg MM is much younger than typical Class I objects, that 
is, it is a Class 0 object or an intermediate phase between Class 0 and I objects (hereafter 
we call it “young protostars”). We should note that the X-ray absorption column is as same 
as those of the Class 0 candidates in the OMC 2/3 cloud (Ah -1—3 xlO 23 cm -2 , Tsuboi 
et al. 2001). 

Xg MM was less X-ray luminous during Obscxo (;$l/30 of ObsxMMi)- Though photon 
statistics are poor, the source looks similarly in activity in another XMM-Newton observation 
on April 2001. This suggests that the X-ray activity was dramatically enhanced during 
ObsxMMi and ObsxMM 2 - Such X-ray variation on long time scales is seen in young stars 
with FU Ori type mass accretion outbursts (Kastner et al. 2004). The X-ray activity during 
ObsxMMi and ObsxMM 2 might relate to a mass accretion episode. 

The slow linear flux increase by a factor of two during —30 ksec in Obs X MM2 is unlike 
magnetically driven X-ray flares in young stellar objects and the Sun, both of which are 
characterized by rapid flux increases within 10 ksec (e.g. Tsuboi et al. 1998, 2000; Stelzer 
et al. 2000; Imanishi et al. 2001). Favata et al. (2003) found a similar event (a factor of 4 
increase during 50 ksec) on the classical T-Tauri star XZ Tau. The event accompanied A H 
decrease as the flux increase so that Favata et al. (2003) interpreted it as eclipsing of the 
emitting region by the accretion stream. We did not find significant variation of spectral 
shape in the hardness ratio of X£ MM (Figure 3). Therefore this model would at least need 
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some modification to apply to Xf MM such as partial covering by an optically thick absorber 
with Nu >10 24 cm -2 . 

Fluorescent iron lines have been rarely observed in stellar X-ray emission. Even a few 
examples during strong flares from Class I sources have EW< 150 eV (Imanishi et al. 2001). 
Large EW of X^ MM again suggests it to be extremely embedded in the molecular cloud 
and/or to have very hot plasma ( kT ~4 keV). Inoue (1985) made a simulation of EW of 
fluorescent iron fine ass umin g possible geometries of surrounding gas (disk -f envelope), a 
power law spectrum of T ~ 1.1 and 1 solar iron abundance for the surrounding matter. 
To apply the result to Xg MM , EW shown in Figure 4 of Inoue (1985) should be reduced 
to ~60% for ObsxMMi and ~80% for ObsxMM 2 because thermal emission of kT ~2.7 keV 
and 4 keV has ~40 and 20% lower X-ray flux around the iron K edge energy (7—9 keV) 
against continuum flux at iron fluorescent line energy (6.4 keV) than a power law r spectrum 
of T ~ 1.1. EW is also proportional to metal abundance. The spectral fitting derives low 
metal abundance of ~0.2 solar, but the result is not likely to apply to the fluorescent line 
because the EW in ObsxMMi exceeds the maximum EW (the case III when X-ray emission 
from a primary source is completely blocked and only fluorescent and scattering emission is 
observed). Circumstellar matter would be reasonably ~1 solar abundance and small metal 
abundance of the plasma should be made by another reason (see discussion in the latter 
paragraph). Assuming solar metal abundance, the best-fit EWs require partial covering of 
the direct X-ray emission by huge covering factor (~60%) for ObsxMMi and small covering 
factor (~3%) for ObsxMM 2 - Intrinsic X-ray luminosity in ObsxMMi is then estimated to be 
log Lx ~31.2 ergs s -1 , which is comparable to L\ in ObsxMM2- 

A consistent explanation for the characteristics described in above two paragraphs is that 
a brightest active spot was behind the proto-stellar core in ObsxMMi and just appeared from 
behind the rim in ObsxMM 2 as a consequence of rotation of the proto-stellar core (Figure 7). 
Certainly, in the case of a solar hot spot, X-ray intensity roughly increases linearly during 
an appearance from the rim (e.g. top panel of Figure 3 of Sterling 1999), and a proto-stellar 
core can be regarded as an optically thick absorber. To be consistent with the observed time 
variation, the active spot should be confined in a small portion on the proto-stellar surface 
(probably within a quadrant) and rotate with a period >2.8 days. Such an X-ray luminous 
small spot cannot be produced by solar type dynamo, and probably produced through mass 
accretion (see Bouvier 1990). The required rotational period is much slower than the break- 
up rotation, speculated for Class 0 protostars from rotational periods of Class I protostars 
(e.g. Montmerle et al. 2000). 

The above picture also suggests that, the intrinsic X-ray activity was constant between 
ObsxMMi and ObsxMM 2 , be. X* MM was in a quiescent phase. The plasma temperature (kT 
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~3— 4 keV) and X-ray luminosity (log L x ~ 31.2 ergs s _1 ) far exceed the typical quiescent 
X-ray activity of Class I protostars, and are comparable to X-ray flares of Class I protostars 
(Imanishi et al. 2001; Shibata & Yokoyama 2002) 5 . This means that youngest protostars can 
continuously produce violent X-ray activity equivalent to flares from Class I protostars. Hot 
plasma temperature requires some acceleration mechanism, perhaps magnetic reconnection, 
in addition to gravitational mass infall. ( Vff ~600 km s -1 of a 1M 0 star can only produce 
plasma with kT ~0.7 keV.) Montmerle et al. (2000) proposed a reconnection of large scale 
magnetosphere by core-disk differential rotation to explain periodic X-ray flares observed on 
the Class I protostar YLW15. A similar mechanism might be applicable to Xg MM though 
magnetic reconnection has to occur continuously to produce constant X-ray emission. The 
model predicts a jet (or outflow) production through the reconnection process (Hayashi 
et al. 1996). Interestingly, Class 0 protostars has an order of magnitude stronger CO outflow 
activity than Class I protostars (Bontemps et al. 1996), which is coincident with the X-ray 
luminosity ratio of Xg MM against typical Class I protostars. 

Our thin-thermal model fit requires a metal abundance of ~0.2 solar, which is mainly 
determined from intensity of the He-like iron line, but the strong fluorescent iron line suggests 
a near solar abundance in the surrounding cold gas. X-rays from pre-main-sequence stars 
show a similar problem regarding small coronal abundance (e.g., Favata et al. 2003). A 
possible solution is that the abundance anomaly localizes around the proto-stellar core. 
A well known phenomenon as such has a relation with first ionization potential (FIP) of 
element, so called the FIP effect. Unfortunately, we can not test the effect on Xg MM because 
abundance of high FIP elements such as argon (0.1< [Z/Zq) At <12 for Model C varying 
argon abundance) or neon (unmeasurable) is not well restricted against iron, a low FIP 
element. Another solution is that the continuum emission includes a non-thermal component, 
as is proposed for p Oph Si (Hamaguchi et al. 2003). 

If the dips in the light curves are really periodic, a possible geometry is that the X-ray 
hot spot confined on the proto-stellar core is blocked by an optically thick gas blob orbiting 
around the core. Orbital radius (r) of the Kepler rotation with rotational period of ~13.9 
ksec is r — 1.2 R Q x ( M/M Q ) J / 3 (f? 0 : solar radius, M: mass of the proto-stellar core). This 
means that a gas blob rotates very close to the stellar surface if it is a solar mass main- 
sequence star. Proto-stellar cores are thought to have several solar radius, and hence the 
derived orbit would be improbable. The dip feature is a riddle. 


5 We should note that Imanishi et al. (2001) use the distance to the p Oph cloud of 165 pc instead of 120 
pc derived from more reliable HIPPARCOS data (Knude & H0g 1998) for a comparison to earlier X-ray 
results. The X-ray luminosity should be divided by a factor of two to compare our result. 
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6.2. What are Xg* 0 and Xw? 

Xg XO has the radio counterpart 10E, but it does not have any IR counterpart (K >19.4 m ) 
though the moderate jV H (~4xl0 22 cm -2 ) is typical of Class I protostars (generally, K <ll m 
in the field). The X-ray and radio emission is steady on short and long timescales (within 
a factor of three change between a day and 2.5 years in X-rays and <15% change in three 
observations between 1985—1998 in the centimeter radio wavelength, Brown 1987; Feigel- 
son et al. 1998; Harju et al. 2001). Its radio spectral energy distribution combined with a 
marginal millimeter detection is consistent with free-free emission (Choi & Tatematsu 2004). 
These results agree with plasma, heated by a collision of a steady jet or outflow with cir- 
cumstellar gas (Tsujimoto et al. 2004). Though X-ray plasma model derived in this paper 
has large uncertainty, log Lx ~29.3— 6 ergs s _1 is comparable to earliest detections of such 
emission in HH2 and L1551 IRS 5 (Pravdo et al. 2001; Favata et al. 2003), while plasma 
temperature >4 keV is much higher than those in HH2 and L1551 IRS 5 (T ^10 6 K) and 
comparable to TKH 8 in the OMC 2/3 cloud (Tsujimoto et al. 2004). Such plasma tempera- 
ture requires an energetic jet with v >1500 km s -1 , which has not been recognized in young 
stars so far. An alternative model with lower kT, higher A H and larger Lx might be more 
reasonable (for example, kT ~1.5 keV, Ah ~ 1.4xl0 23 cm -2 and log Lx ~30.3 ergs s -1 , 
which are acceptable to the Xg XO spectrum though it still needs v >900 km s -1 and L x 
comparable to stellar origin). On the other hand, the X-ray emission does not show any 
spatially extended feature, suggesting that the X-ray emitting (perhaps bow shock) region 
should be confined smaller than 1" (~170 AU). If this hypothesis is correct, Xg MM (~340AU 
projected distance) is the closest and therefore the most plausible jet driving source. The jet 
would flow to the south-east direction, probably aligned with the rotational axis of Xg MM 
(Figure 7). We should note that Class 0 protostars drive larger momentum CO outflows 
than Class I protostars (Bontemps et al. 1996), and therefore would produce hotter plasma 
at the jet colliding region. 

Xw has a bright near-IR counterpart IRS 7 as well as the centimeter radio counterpart 
10W, and therefore might be a stellar origin. The Ah ~2— 6xl0 23 cm -2 is comparable to 
%xmm, w j 1 j] e Riband magnitude (~12.2 m ) is comparable to Class I protostars in the 
field. X w (IRS 7) might be a Class I source with a large inclination angle. On the other 
hand, the ASCA hard X-ray source exhibiting a strong X-ray flare in 1994 (Koyama et al. 
1996) has been suspected to be IRS 7, but IRS 7 (Xw) has an order of magnitude larger Ah 
than the ASCA source (Ah ~4.3 x 10 22 c,m~ 2 ), and hence would not be the counterpart. 
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7. Summary &c Conclusion 

We discovered three extremely embedded X-ray sources in the IRS 7 core in the XMM- 
Newton and Chandra observations. In particular, Xf MM has not been recognized at any 
wavelengths, which means that X-ray observations can be a powerful tool to search for 
young protostars. Thanks to its vicinity to the Sun (d ~ 170 pc), the large effective area of 
XMM-Newton and a rare opportunity to catch an active state, we could obtain around ~2,000 
photons from the young protostar X£ MM , which is around 40 times better than those obtained 
from X-ray emitting Class 0 candidates so far. The upperlimit to the R-band luminosity 
was restricted to — 1/150 of the OMC 2/3 protostars (d ~450 pc and completeness limit 
~16 m , Tsujimoto et'al. 2003). Xg MM would thus represent the most reliable X-ray detection 
of an example for the youngest protostars, and the best characterized in both IR and X-ray 
properties so far. The multiple X-ray observations and excellent photon statistics obtained 
with XMM-Newton enabled us to study the X-ray emission mechanism in further detail. 
The X-ray activity might be driven by strong mass accretion in the earliest phase of stellar 
evolution. The activity scale (fcT~3-4 keV, log L x ~31 ergs s _1 ) was actually far stronger 
than constant X-ray activity of Class I protostars, which would probably connect to the 
strong CO outflow phenomena seen around Class 0 protostars. The apparent time variability 
observed during the XMM-Newton observations could be produced by slow rotation of the 
proto-stellar core with a period of >2.8 days. The period is much slower than the break-up 
rotation, speculated for Class 0 protostars from rotational periods of Class I protostars (e.g. 
Montmerle et al. 2000). This may mean that angular momentum is effectively removed from 
accreting gas during the mass accretion process. The weak metal abundance feature of ~0.2 
solar was also seen in young protostars. This suggests that local abundance anomaly around 
the proto-stellar core can develop within ~10 4 years or the feature is apparently made by 
fundamental X-ray emission mechanism such as non-thermal process. On the other hand, 
Xe XO has some characteristics of a jet induced plasma. Xg MM and Xg XO might form an 
X-ray jet driving and induced system. 

The evolutionary status of Xg MM as whether a Class 0 protostar or an intermediate 
phase between Class 0 and I protostars should be elucidated by constructing the IR and 
radio SED. For the purpose, we need high spatial resolution and sensitive observations in 
IR and radio wavelengths. Ground based 10m class , Spitzer and ALMA telescopes will be 
good instruments to reveal the nature of this mysterious X-ray protostar. Detection of a 
jet feature and measurement of the jet direction would help to know the inclination of the 
rotational axis as well. On the other hand, revisit of the field with X-ray observatories with 
a long exposure could detect periodical X-ray time variation of Xg MM by the core rotation, 
restrict plasma nature of Xg XO and Xw with better photon statistics, and discover other 
X-ray emitting young protostars, which would be hidden inside gas clumps in the field. 
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Table 1. Observation Logs 


Observation 

Observatory 

Seq. ID 

Date 

Exposure 

(ksec) 

Aaxis 

Shift {Aa.AS) 

ObsxMMi 

XMM-Newton 

146390101 

2003 Mar. 28 

14.9/21.5 

6' 

(0." 4, -1." 0) 

ObSxMM2 

XMM-Newton 

146390201 

2003 Mar. 29 

18.1/24.0 

6' 

(0." 2, -1." 3) 

Obscxo 

Chandra 

200017 

2000 Oct. 7 

19.7 

2' 

(0." 03, 0." 23) 


Note. — Exposure: pn/MOS for XMM-Newton. A axis'- Off-axis angle of the IRS7 region. Objects 
used for the position correction: IRS2, IRS5, HBC677, CrAl, ISO-CrA136, HH101 IRS 4, ISO-CrA137, 
TY CrA and HD 176386. IRS1 and ISO-CrA134 were also used for the position correction of Obscxo- 


Table 2. Detected Sources 


Source 

Abbr. 

Observation 

( a 2000 * $200o) 
(h m s, d ' ") 

Error* 

Net cnts 
(cnts) 

K-band 

(mag) 


Counterpart 

XMMU J1901 56.3-365726 

vXMM 

x E 

ObsxMMi 

19 1 56.29, -36 57 26.5 

1." 7 

467.7 

>194 





Obs X MM 2 

19 1 56.32, -36 57 26.4 

1." 2 

1611.3 




XMMU J 1901 55. 3-365721 

X^ MM , Xw 

° bs XMMlfc2 

19 1 55-26, -36 57 20.8 

-w3" 

U4.0 b 

12. 2 C 

10W d 

, Source 4 e , IRS7 

CXOU J1901S6.4-365728 

xg xo 

Obscxo 

19 1 56.43, -36 57 27.7 

1." 0 

19.6 

>19.4 

10E d 


CXOU J 1901 55. 3-365722 

Xff®, Xw 

Obscxo 

19 1 55.34, -36 57 21.6 

1." 0 

12.7 

12. 2 C 

10W d 

, Source 4 e , IRS7 


Note. — Net cnts: EPIC pn 4- MOS for XMM-Newton. 

*90% confidence error. 
b Net counts in ObsxMMi- 

c Reference of the K x magnitude of Rl (IRS7) to Wilking et al. (1997). 
d Reference to Brown (1987) and Feigelson et al. (1998). 

** Reference to Choi ic Tatematsu (2004). 


Table 3. Fitting Result of the Xg MM Light Curves 


Observation 

Binning 

(sec) 

Model 

Constant 
(10 -2 cnts s _1 ) 

Linear 

(10~ 2 cnts s _1 day -1 ) 

X 2 /d.o.f. (d.o.f.) 

ObsxMMi 


con. 

0.93 


1.63 (13) 

ObSxMM2 


con. 

5.3 


1.96 (26) 



con.+lin. 

3.7 s 

9.3 

1.21 (25) 


a Count rate at TJD = 12727.8564 day. 


Note. — con.: constant, lin.: linear 










Table 4. Fitting Results of the Spectra 


- 20 - 


5 ~T' 
L 


Eh > 


a ; 

^ C4 


(N CO CS ( 


X ifi 
if lO 

I I 


if if oo if 
cb rb cb cb 


O .-I r-< *f tO cs n o 

Cob’T'TMN’Tb 
lO lO LO U) A m lO lO 

l I I I I I I I 

O O k> O co oo oo © 

cb if cb o* -i cb cb 
to © to © to © uo to 

n h * h io .-< n oo 

cb if cb if cb cb cb cb 

to to to to to to to to 


boo 

'll ~ 

T T „ x 
oo.q h Jib 


LO 04 > > 

6 o 6 II II 


CO co 

o d 


0 . • - 

-* to v to ^-v 

1 I I I I 

M 'f CO O) 

oi cb ci o co 

« i; n o q ^ 

© if cb if cb A 


. — « co oo 
” oo cb 
oo ci co 

7 1 1 

I <o m 
«? rb cb 
o cs cs 


7 i 


h h H H Eh h 

H H N i-H .-H r-l 

< O O QHfc 


JO JO -O 

O O O 


xa 

X 


O 0 o 


0 I 0 

okxSoS 

XXX 


£ & 


3 2 

.5 tj 



- 21 - 


R.A. ( J2000 ) 

10s 1 9h02m00s 50s 40s 19h01m30s 




TV CrA 

m 





m 

HI) 1763S6 


K CrA 






IRS? 



IRS" 

region 

IKS1 


iks: 

ISO CrA 134 

\ , 



HBC 677 


/ 

ISO Cr A 136 


-36:52:00 


54:00 


56:00 


58:00 


a 

P3 

n 


© 

© 


-37:00:00 


02:00 


Fig. i _ XMM-Newton true color image of the R CrA star forming region (Obs X MMi + 
0bsxMM2)- The image is color coded to represent hard band (3—9 ke\ ) to biue, medium 
band (1-3 keV) to green and soft band (0.2-1 keV) to red. The dotted square shows the 
IRS 7 region (field of view of the top panels of Figure 2). Class I protostar: IRS 1, IRS 2 and 
IRS 5. Herbig Ae/Be star: R CrA. TY CrA. HD 176386. Weak-lined T-Tauri star: CrA 1. 
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Fig. 2. — Magnified images of the IRS 7 region. The top left panel shows the 3—9 keV EPIC 
pn + MOS image of XMM-Newton (ObsxMMi + 0bsxMM2)> the top right panel shows the 
0.5—9 keV image of Chandra (Obscxo) and the bottom right panel shows source positions 
with their error circles including sources in near-IR and radio wavelengths. Diamonds show 
XMM-Newton sources, crosses Chandra sources, squares centimeter radio sources (B number, 
Brown 1987, Source positions are referred to Feigelson et al. 1998), x millimeter radio sources 
(GT number, Choi & Tatematsu 2004), a triangle near-IR source (IRS 7). Error circles are 
shown in dot-bar line for X-ray sources and dotted line for radio sources and solid line for 
near-IR source. Error circles of X-ray sources for the 2MASS frame is smaller than the above, 
and those of Xg MM (0bsxMM2) and Xj? xo do not overlap with each other. 
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Fig. 3. — Light curves of X^ MM in the 2—10 keV band ( top panel) and the hardness ratio 
defined as ” count rates in the 5—10 keV band / count rates in the 2—5 keV band” ( bottom 
panel). Points with error bars in the top and bottom panel are made from summed data sets 
of pn plus MOS, while grey diamonds and an error bar at 7.85 day in Obs X MM 2 show fight 
curves of MOS (1+2). The horizontal axis is truncated Julian day (TJD) — 12720. The 
black arrows show timing of the dip feature and the grey arrow does a possible dip timing 
assuming the period of 13.9 ksec. Bins in the top panel have 2 ksec (Obs X MMi) and 1 ksec 
(Obs XMM2 ), and bins in the bottom panel have 4 ksec (Obs X MMi) and 2 ksec (Obs XM M2)- 
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Fig. 4. — EPIC pn (left) and MOSl+2 (right) spectra of Xg MM in Obs X MMi (black) and 
Obs X MM2 (red). Solid lines show the best-fit model of the simultaneous fitting with EPIC 
pn and MOS. Dotted lines show the constant (soft) component, barred line the variable 
components, and dot-bar lines the Gaussian component. The bottom panel shows residual 
of the x 2 fit- 








Fig. 5. — Chandra spectra of Xg XO (left) and X w (center), and XMM-Newton spectra of X 
in ObsxMMi (right, black-pn, red-MOS). Solid lines show the best-fit models (Model D, 
and E in Table 4 from left). 
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Fig. 6. — The TF-band image of the IRS7 region obtained in 1998 ( top left : original image, 
bottom: fit of the original image by the median-filtering method, top right: subtracted image). 
X| MM are at the cross of two long (horizontal and vertical) bars. Xg XO (left) and X^ 0 (Xw, 
right ) are at the crosses of short bars. 
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Fig. 7. — A possible geometry of the proto-stellar core, jet and observer . The wavy lines 
show X-ray emission (solid black: direct X-rays, gray: fluorescent and scattering X-rays, 
dotted black: direct X-rays are partially covered). 


